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I. INTRODUCTION 


The output signal of a radiometer can be expressed in terms 
of the distribution of incident spectral radiance, the instrument 
parameters and relative geometry between the source of the radiation 
and the instrument. The general formulation of this relation is 
commonly referred to as measurement equation [^ij. The application 
of this measurement equation to the solution of practical measurement 
problems usually requires some simplifications to obtain a tractable 
computational model of the measurement process. The primary objectives 
of this study are: (1) to derive (simplified) formulations of the 

measurement equation suitable for evaluating effects of the spatial 
(or angular) response of narrow f ield-of-view (NFOV) scanning radio- 
meters; and (2) to evaluate effects of the measurement process that 
are inherent in the reconstruction of a continuous radiance field 
from discrete radiometric measurements. 

In particular, it is shown that appropriate approximations 
to the measurement equation lead to a convolution integral of the 
spatial distribution of the spectral radiance and of the radiometer 
response, and that the accuracy of continuous recons tructions of 
radiance fields from discrete radiometric measurements is subject 
to degradation due to blurring and aliasing. It is proved (in the 
Appendix) that the aliasing errors that are generated if spatial 
details have been undersampled, can be treated as statistically 



independent noise if the radiance field is a (real) Gaussian process. 
Degradations due to aliasing and blurring are evaluated for a random 
radiance field with a Wiener spectrum that is representative of a 
wide range of scenes, and for spatial responses that are typical 
of televisions and radiometers. In addition, this study also 
evaluates the information capacity representative of these spatial 
responses . 

This study is part of an undertaking to develop comprehensive 
analytical and computational measurement models for the Earth 
Radiation Budget Experiment (ERBE) . The results obtained in this 
study have been used in the ERBE NFOV design and performance trade- 
off studies. 

The measurement of the Earth's radiation balance to refine 
previous estimates has attracted continuing attention over the 
last two decades. These efforts, which included experiments on 
Explorer VII and the Nimbus satellites, have led to the Earth 
Radiation Budget Experiment sponsored by the National Aeronautics 
and Space Administration (NASA). The ERBE instruments include 
wide and medium field-of-view (WFOV and MFOV) radiometers, as well 
as narrow field-of-view (NFOV) scanners. While the wide field-of- 
view instruments are useful in obtaining global averages, the 
narrow field-of-view instrument provides local measurements of 
radiation, which are useful in the accurate determination of spatial 
(zonal and regional) variations of the radiant exitance at the top 
of the atmosphere, and can be used in the estimation of the directional 
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models of the radiance field. To be useful, an accuracy of 1% for 

the global mean albedo and radiant exitance (a and resp.), and 
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an accuracy of ±10 W/m for zonal and regional variations in the 
short wave (reflected) and long wave (emitted) radiant exitance 
are generally required for climate monitoring, with more stringent 
accuracies required for climate theory applications. Thus, it is 
necessary to determine the accuracy of the measurements and investigate 
the trade-offs involved between the design parameters and the errors 
introduced in the measurements of the radiant exitance at the top 
of the atmosphere. 

This work considers the analysis of aliasing error and blurring 
(or the loss of spatial resolution) introduced into a reconstruction 
of the radiant exitance when the NFOV scanner output is undersampled. 
Through mathematical modeling and quantitative analysis, it describes 
the trade-offs between design parameters such as the sampling 
interval, detector aperture shape, continuous-scan versus step-scan, 
etc. , and the amount of aliasing error and blurring introduced into 
a reconstruction of the radiant exitance from the sampled digital 
radiometer output. The effects of smoothing the estimates, and of 
varying levels of high spatial frequency content in the radiant 
exitance at the top of the atmosphere are also considered. To 
obtain results useful in a more general context, when possible, 
normalized parameters were used. The degradation of image quality 
obtained with line-scan devices such as optical-mechanical scanners 
and television cameras has received considerable attention [^2^ - 
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jj.Oj . In these investigations, the emphasis has been placed mainly 
on the detection, recognition and reproduction of specific detail 
and objects. The work presented here considers a detailed para- 
metric study of the expected magnitude of aliasing errors contained 
in the reconstruction of random radiance fluctuations typical of 
natural scenes, which has been considered only in illustrative 
examples in the literature. 

The block diagram in Fig. 2 shows the main components of the 
NFOV scanner considered. The radiant exitance M(x,y) at the top 
of the atmosphere passes through a lens and heats the detector, 
thus altering its temperature. The change in the temperature is 
sensed and converted to an electrical signal which is amplified 
and passed through an electrical filter to reduce the noise due 
to the detector and aliasing error, before being sampled, digitized 
and transmitted to earth by telemetry. The sampled signal is then 
reconstructed to obtain a continuous signal. In most cases, due 
to limitations in the telemetry transmission rate and the frequency 
content of the signals, the reconstructed signal differs from the 
original one; this is referred to as insufficient sampling. The 
error in the reconstructed signal is called aliasing noise or 
error. In general, the smaller the sampling interval, the less 
the aliasing error introduced. However, as can be seen from the 
results (see section III) the choice of design parameters such as 
the shape of the detector field of view, the electrical filter, 
the use of a stepping or a continuous scan system, etc., is of 
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considerable importance in reducing the aliasing error. 

In section II, the relative geometry of the orbit, the top 
of the atmosphere and scanning pattern (see Fig. 1) are used to 
show that the output of the detector can be modeled or approximated 
by passing the radiant exitance through a space- invariant system. 
Section IIA describes analytical models of the aliasing error and 
its dependence on the system parameters. Section IIB considers 
the information density transmitted through the system as a design 
criterion. Section III describes the results of the parametric 
study to show the trade-offs involved quantitatively. 


5 



II. NFOV SCANNER MODEL DEVELOPMENT 


Consider a satellite on a circular orbit about the earth, with 
radius r + h, where r is the radius of the sphere representing the 
top of the atmosphere, and h is the altitude of the satellite from 
the top of the atmosphere, as shown in Fig. 1. The subsatellite 
points of the orbit form a great circle on the top of the atmo- 
sphere. Considering this circle as the equator, a coordinate 
system with components of longitude (cj)) and colatitude (0) can be 
defined, to represent the points at the top of the atmosphere. 

Thus, the equator and constant colatitude curves represent the 
along- track direction, while the constant longitude curves re- 
present the scan direction. 

Suppose that, at time t, the subsatellite point has longitude 
(f)^, and that the scan point (the intersection of the optical axis 
and the top of the atmosphere) has colatitude 0^, then the detector 
output, m(0^, radiation from the top of the atmosphere 

can be modeled as (e,g.,. see Chapter 9) 

m(e, , ij) ) = A // /“ L,(6, (fi, 6’, (}>’)S,(X, <1^) dXdA' (1) 

1 1 ^ A p2 

where A is the area of the detector, (6 » > 9’» 4^') is the 

spectral radiance at the top of the atmosphere at the point (0, 0) 
in the direction (0’, (p'), is the sensor response for 

radiation arriving from the direction (x» of wavelength A, p is 
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the distance between the point (0, 4>) and the detector, and is 
the ’’footprint" at the top of the atmosphere for the scan point 
(01, 4*2^)* This is based on assumptions that the instrument is 
operating in a stable and linear region of its operating range, 
that the effects of radiation arriving from other sources have 
been accounted for by ground and in-flight calibration procedures, 
that the sensor time constants are small compared to variations in 
the spectral radiance, and that the linear dimensions of the detec- 
tor surface are much smaller than the distance p. 

For a narrow field of view sensor, Eq. 1 can be approximated by 
a convolution as will be shown in the following; this allows the use 
of Fourier theory in the analysis of the measurement process of the 
NFOV scanner. The scan point with longitude (j?!, colatitude 0i, and 
radius r can be represented by a three component vector, 
the rectangular coordinate system shown in Fig. 1, as 


Z 2 = r 


sin0i cos(|)i 
sin©! sin4>i 
cos0^ 


( 2 ) 


The position of the detector, say z^ can be expressed as 


z^ = (r + h) 


COS( 


sincf). 


( 3 ) 


Now let z be a point in the footprint A| about the scan point 


7 


(©I* ^1 ” colatitude and longitude of 

z, so that 


z = 


r 


sin(0^ 
sin(0^ 
cos (0j^ 


a) cos((J)^ - 3) 
a) sin((j)^ - 3) , 
a) 


For a narrow field of view instrument, the footprint consists of 
a small area about the scan point, so that the values of a and 3 
are small. Thus, using small angle approximations. 


z - 


r 


(sin0^ - a cos0^) (coscf)^ 
(sin0j^ - a cos0^) (sincj)^ 
cos0^ + a sin0^ 


+ 3 sincf)^) 

- 3 coscj)^) . 


Neglecting the second order terms in (5), 


z - z„ -(ra u„ + r3 sinO. u,), 
2 6 1 (p 



COS0^ COs4)^ 


- sin<f)^ 

“0 “ 

cos0^ sincf)^ 

• > = 

coscf)^ 


- sin0- 


0 


1 _ 




It should be noted that u^ and u^ are orthogonal unit vectors which 




lie in the plane tangent at the scan point, z^. Thus, from (6) it 
is seen that the points in lie (approximately) on the tangent 
plane. To illustrate the approach, consider a circular field of 
view of 6 radians, then 



^ < X £ <5} = {(0^-a, 4)^-3) 1 0 < sinx < sin6} (8) 

For a given scan point, the angle X is a function of a and 3; 
thus, the constraint on X (S) imposes a constraint on a and 3» 
which defines the footprint more explicitly. To determine the 
latter constraint, first note that X is the angle subtended by the 
optical axis and the line connecting the detector to an arbitrary 
point in the footprint; i.e., x is the angle between the vectors 
^2 ~ z - 2^. Hence, 


sin X = 


(z^-z^) X (z-z^) 


II V^l It 


z-z. 


(9) 


where |j j] denotes Euclidian vector norm, and x denotes the vector 
product (cross product) . Noting that 


z-z^ = (z-z^) + (z^-z^) , (z^-z^) X (z^-z^) = 0 


( 10 ) 


(9) can be expressed as 


IKvV ’ 


( 11 ) 


where 


P = II II » P2 


^2 ^ 1 II * 


( 12 ) 


Using the definition of the vector product. 


sinx = 




P2IIZ-Z2 


T 


( 13 ) 


where the superscript "T" denotes the transpose. From (6), we 
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obtain 


2 -Z 2 II ||aug + 3 sin0^ ||^ 


2 2 2 2 
= r (a + 3 sin e^) , 


Manipulating (2) and (3), 


+ (r+h)^ - -2r(r+h) sinG^, 


By direct computation, it can be found that 




(z^-Zj^) Uq = “(r+h) COS0^ 


T T 

(z^-z^) (z-z^) = -rCz^-Zj^) (a Uq + 3 sinG^ u^) 


= - ra(z -z ) u„ = r(r+h)a cosG- , 
z i y 1 


Substituting (15) and (20) into (13), and collecting terms, 


sin X ~ 




2 2 2 2 2 
(r+h)^ cos G^) a + P 2 sin G 


P P 


(r+h) sinG^)^ + p^ sin^G 


.=1 


(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 


( 21 ) 

( 22 ) 


where (16) has been used to obtain (22). A further simplification 
can be obtained by noting that the angle, G|, subtended by z^-z^ 
and the normal to the tangent plane at (0^^, (f)^) , satisfies 
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cos 0 


"2 


2 j_-Z 2 11 r p. 


— ((r+h) sin 0_ - r) . 

p2 1 


f T. 2. 

(Z2 z^-r ) 


(23) 

(2A) 


Substituting (24) into (22), and simplifying 


sin X 


2 2 2 2 
cos 0j^ a + sin 0^ 3 


Since the footprint, A^, corresponds to a small area about 
the scan point, or alternately since a and 3 are small for points 
in the distance p can be approximated by p^. Thus, the foot- 

print can be described by 

2 2 2 


{(0^-a, (j)^-3) 1 005^^0^ + sin^0^ 3^ < -^ 5^} 

r 


(25) 


(26) 


Thus, for a narrow circular field of view, the footprint 
corresponds to a nearly elliptical surface on top of the atmosphere. 
At nadir, the footprint becomes circular as a limiting case, since 
0^ = 7t/2 and = 0. 

It is sometimes more convenient to use distance related 
variables than latitude and longitude. Thus, consider the following 
change of variables 


X = r({), y = r0, x'= r3, y’= ra. 


( 27 ) 


In terms of these variables, (1) can be rewritten as 


11 



m 


(x,,y,) = A // /” L,(x, y, 0', <()’)S,(x, « dXdA' . 

Ai 0 A A Z 


Now let 


X = x_ - x’ , y = Yt - y’ ; 


Then dA’ is given by 


dA’ = sin0 d6 d({) = sin0 dx’dy' 


Define the indicator function of A| as 


Ij^Cx’ ,y' ) = 


. . 2 . ,2 , 2 „, ,2 ^ 2 .2 

1, sin 6^ X + cos y ^ P 2 ^ 


0, otherwise 


Substituting (29), (30) and (31) into (28), 

00 oo 00 I 

m(x^,y^) = A /^ /^ / L^(x^-x’, y^-y’, 0’, 4>’ ) |^S^(X,4^)I^(x' ,y 


COS0' sin0 j 1 j I 
r dA dx dy 


Over the footprint A|, 


sin0 = sin(0^-a) = cosct - cos0^ sina - sin6^, 


0' 0|, r - 4>;, P - p2- 


With these approximations. 


A COS0' sin0 oo c» co 

mCxj^.yj^) = 2 L iZ { ’ yi~y'’ ®’> 


js^(X,'J')Ij(x',y') cosxjdXdx'dy’ 


(28) 

(29) 

(30) 

(31) 

’ ) COS 

(32) 

(33) 

(34) 

(35) 
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To obtain (35) as a convolution, it is necessary to separate 
the wavelength, X, from the other variables. Thus, we assume that 
the spectral radiance at the top of the atmosphere admits a spectral- 
directional model of the form 

L;^(x, y, 6’, (f)’) = ^ M(x,y) R^(0 ’ , (j)'), (36) 

where M(x,y) is the radiant exitance at the point (x,y) at the 
top of the atmosphere. It may be noted that the desired form 
(convolution) can be obtained with models somewhat more generalized 
than the one shown in (36); however, these generalizations will 
not be pursued here. We further assume that the spectral and 
angular characteristics of the sensor are independent of each other; 
i.e. , 


S^(X,^) = S(Xy^) 


(37) 


Rearranging (25) in terms of x' , y’ 


X(x’ , y ' ) - sin ^ — {sin^G x’ ^ + cos^G' y 
p/-, 1 1 


(38) 


Similarly, it can be shown that the angle i{j is also a function of 
x' and y'. Substituting the expressions for x ^nd ip, we define 


5^(x’, y') = S(x(x’,y'), i}j(x\y')) I^(x’,y’) cosx(x’,y') 


(39) 


Substituting (36), (37) and (39) into (35), 
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A COS0’ sin0- oo oo c» i 

m(x^,yj^)= C L. f — M(x^-x’.yj^-y’)Rj^(e|,(t.ps^ 


-00 -0O Q IT 


S^(x’,y*) dXdx'dy’ 


00 oo 

m(Xi,yi) = M(x^-x’, y^-y ' )S^(x’ ,y ’ )dx' dy ’ 


A COS0' sin0. 

K 2 r 

TTP 2 

Thus, it is seen that m(x^,y^) can be obtained (approximately) 
by a convolution as shown in (41). However, and still depend 
on the scan point (x^,y^^). Thus, consider a "patch" about the 
point (XqjYq) such that for any scan point (x^,y^) on the patch 


sxne^ “ s“9o’ ®i = ®0’ ‘*'1 = ‘^0’ ^2 “ ^0 


Then, for any point the patch 

A COS0' sin6 

^0 r— ^ r 

xPq 

OO JX> 

m(x^,y^)=K^ M(x^-x' , y^-y') S^(x’ ,y’ )dx’dy ’ 

Thus, (45) is the desired form for the measurements, showing that 
a narrow f ield-of-view scanner can be interpreted as a linear, 
spatially invariant system where SQ(x’,y') is the point spread 
function of the system. It is of interest to note that the point 
spread function is composed of two parts. The effect of the field 
of view is incorporated by lQ(x’,y’) which vanishes except for 


(40) 

(41) 

(42) 

(43) 

(44) 

(45) 
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points within the field of view. It should be noted that the effect 
of the partially obstructed field of view can be included by appro- 
priate choice of Iq(x’ , y’). The angular characteristics of the 
detector are incorporated by the product S(x, cosx* Ideally, 
S(x>^) would be constant, at least over the field of view; and, for 
a narrow f ield-of-view sensor, cosx - 1. Thus, for an ideal 
sensor, the point spread function would be consent for points 
within the field of view and zero for other points. 

It should be noted that, in terms of the coordinates x' , y' , 
the point spread function S^(x',y’) varies according to the patch 
considered. For the circular field of view considered, it can be 
seen from (31) that the footprint will be circular at nadir, but 
will become elliptical at other points, depending on the colatitude 
of the patch considered. Furthermore, note that the conditions 
required for a patch (Eq. (43)) are independent of the longitude 
or equivalently x^ , but depend on the colatitude 6^ . Thus, a 
patch could be defined as an annular surface between two equal 
latitude curves. 

A further point of interest is the time dependence of the 
radiant exitance, say Mq(x, y, t) , at the top of the atmosphere. 
Suppose that at time t, the scan point is given by (x(t), y(t)), 
and define 

M(x(t), y(t)) AM^(x(t), y(t), t) . 

If the change in the radiant exitance during the time required for 
the satellite to move a distance corresponding to its field of 
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( 46 ) 



view is negligible, then (45) remains unchanged when M is defined 
by (46). 

m(x^(t), = Kg M(x^(t)-x' yj^(t)-y') Sg(x’,y') dx'dy'. 

Thus, the convolution remains unchanged, and the time dependence 
of the signal m can now be seen explicitly in (47) . 


A. Aliasing Error Model 


The block diagram in Fig. 2 shows the main components of the 

ERBE NFOV scanner which is then followed by a reconstruction of the 
sampled signal. Thus, the radiant exitance M(x,y) at the top of the 

atmosphere passes through a lens and heats the detector, thus 
altering its temperature. The change in the temperature is sensed 
and converted to an electrical signal which is amplified and 
passed through an electrical filter to reduce the noise due to the 
detector and aliasing error, before being sampled, digitized and 
transmitted to earth by telemetry. The sampled signal is then 
reconstructed to obtain a continuous signal. In most cases, due 
to limitations in the telemetry transmission rate and the frequency 
content of the signals, the reconstructed signal differs from the 
original one; this is referred to as insufficient sampling. The 
error in the reconstructed signal is called aliasing noise or 


( 47 ) 
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error. In general, the smaller the sampling interval, the less the 
aliasing error introduced. However, as can be seen from the results 
(see section III) the choice of design parameters such as the shape 
of the detector field of view, the electrical filter, the use of a 
stepping or a continuous scan system, etc., is of considerable 
importance in reducing the aliasing error. Furthermore, the re- 
duction in aliasing error obtained as the sampling interval is 
reduced is not uniform; so that in some instances a small reduction 
in the sampling interval can reduce the aliasing error noticeably 
more at some points than at others. Thus, in order to make the 
trade-off between design considerations (cost, complexity, etc.) 
and aliasing error, it is necessary to make a parametric quantitative 
study of the relationships between the various parameters involved. 

It has been shown (Eq. (45)) that the effect of the detector 
can be expressed by a linear space- invariant system, so that the 
output can be described as a convolution of the input with the 
point spread function of the system. To establish the notation 
which will be used, let g(x,y) be an integrable function; then 
its Fourier transform g(v,o)) is given by 


oo oo 

g(v,w) = g(x,y) 

OO oo 

g(x,y) - g(v,o)) 


-i27T (Vx+ioy) , . 

e ■' dxdy. 

i2TT(vxfojy) , - 

e ■' dvdo). 


(48) 

(49) 


In these expressions, g(x,y) can also be interpreted as a second 
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order stationary (i.e., homogeneous and isotropic) random process; 
for a more detailed treatment of these definitions, see Appendix A. 

It is well-known that the Fourier transform is linear and that the 
transform of a convolution is the product of the transforms of the 
convolution components [l2] , [l3] . Thus, the transform of the 
detector output given by (45) is given by 

m(V,03) = M(V,w) (50) 

where M(V,w) is the transform of the radiant exitance at the top of 

atmosphere and T^(V^a)) the transform of the point spread function 

T,(x,y) defined as 
d 


Td(x,y) = Kq S^(x,y) ; 


(51) 


we will refer to t^(v^ 0 }) as the frequency response of the detector. 

The effect of the lens has been neglected in (50) , as the diffraction- 

limited spatial frequency response of the lens does not usually 

affect the frequency response of relatively low-resolution radiometers. 

Simulations, including the lens frequency response, during the early 

stages of this study, also led to the above conclusion. Thus, the 

effects of the lens will be neglected in the following. 

In a step-scan system, the output of the detector is effectively 

sampled by the stepping process; i.e., if the scan points are 

(x.,y,), then the output of the radiometer at these points is 

m(x.,y ). Thus, after a period is allowed for transient effects to 
J k 

die out, the detector output can be sampled, digitized and transmitted 
to ground . 
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When a continuous scan system is used, the detector output is 

amplified, then passed through a low-pass electrical filter before 

it is sampled. The low-pass filter serves two purposes: to reduce 

the amount of noise, and to lower the high frequency content of the 

signal before it is sampled, thus reducing aliasing error. Hence, 

a continuous-scan system provides a method of reducing both noise 

and aliasing error which is not available in step-scan systems. 

The amplifier gain will be assumed to be unity, for simplicity. 

Now, if the scanning speed is much faster than the along- track 

speed, then the time rate of change of the detector output follows 

the variations along the scan direction. In fact, let h (t') be 

e 

the impulse response of the low-pass filter; then the output of the 
filter can be expressed as 


s(x(t) , y(t) ) = 

h (f) 
e 

m(x(t- 

t’), y(t-t 

’)) dt' 

( 52 ) 

OO 

== Zoo 

h^(t') 

m(x(t) 

, y(t-t’)) 

dt' 

( 53 ) 

OO 

= Zoo 

hjt') 

m(x ( t) 

, yt-yt’) 

dt’ 

( 54 ) 

II 

T^(y-) 

m(x(t) 

, y(t)-y’) 

dy’ 

( 55 ) 

where y is the scanning 

speed , 

and 





y(t) = yt , y’ = yt' , ( 56 ) 

T^(y’) = h^(y'/y)/y. 


( 57 ) 



Thus, the effect of the electrical filter can be approximated by a 
convolution in the scan direction. Taking the Fourier transform 
of s(x,y), we obtain 

s(v, (jo) = T, (v,(jo) t (oo) M(v,co) = t(v,(jo) M(v,w) , (58) 

Q e 

✓v 

where t(V,(jo) is the combined frequency response or MTF of the radio- 
meter till the filter output. Note that the effect of the lens 
results in the MTF 


T(V,M) = Tj^(V,0o) T^(v,u) T^(lO). 


(59) 


The output of the filter is then sampled uniformly in time, 

corresponding to sampling intervals X and Y along the x and y 

directions, respectively. This results in samples of the form 

s(x,,y, ); we denote the sampled signal by s(x,y; X,Y) . The Fourier 
3 k 

transform of the sampled signal is then given by 


s(v,lo; X,Y) 




(60) 


(61) 


If the Wiener spectrum (or power spectral density) of the radiant 

-A. 

exitance is denoted by oj) , then the Wiener spectra cf) (V, oi) and 

$^(v, co; X,Y) of s(x,y) and s(x,y; X,Y) , resp. , are given by 


(f)^(\;, CD) = |t(V,(D)| 0j^(V,(D') 


(j)^(V,(D 




j =-oo k=-co 


_ Jl 


^sV"^ - ^ 


(62) 

(63) 
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It is seen that the spectra of the sampled signal and the 
continuous signal differ by the aliasing error spectrum <p (v,cxj; X,Y) 
Under sufficient sampling conditions, the aliasing term vanishes, 
and it is possible to reconstruct the continuous signal from its 
samples with no error. When insufficient sampling occurs, this 
term is non-zero, and any reconstruction contains errors , [^18^ • 

Consider the reconstruction 


OO C50 


r (x 


»y) = y J y J sin c-^^^ sin c^"^^ 


j =-“ k=-«= 

which corresponds to the reconstruction filter with frequency 
response 

1 , Ixa^I < |yoj| < 

0 , otherwise 


T (-V,(jj) = n(XV, YO)) = 
R 


From the sampling Theorem [18] , it follows that 


r(v, oj) = Tj^(v, w) s(\), w; X,Y) 

= s(v,[jq) + Tj^(v,co) s^(v,03; X,Y) 

= r (V,w) + r (v,oj; X,Y) . 
s a 


(65) 

( 66 ) 

(67) 

( 68 ) 
(69) 


Taking the inverse Fourier transform of (69) 



It is seen that the reconstruction contains a part that is due 
directly to the signal, and a part that is not desired, called the 
aliasing error. In fact, it can be shown that if M(x,y) is a (real) 
gaussian process, then the signal and aliasing parts of the recon- 
struction are statistically independent and stationary processes 
(see Appendix A, Th. 3). Thus, the aliasing error term can be 
treated as additive noise. 

If the useful frequency band of the signal before sampling 

s(x,y) is considered to be B = {(v, 03 )l |Xv[ < ht lYoij < h) t then, 

with T (v,oo) given by (66), r (v,a)) is precisely the useful part of 
R s 

the signal. Now, since r (x,y) and r (x,y) are stationary processes, 

S Si 

their variances do not change with the position (x,y) . 


9 2 00 <0 

O = E(r (x,y)) = / / 


Tj^(v,0)) 


0g(v>oj) dvdoj 


(71) 


.1/2X A/2Y 
-1/2X -1/2Y 


$^(v,co) 


dvdco 


(72) 


2 2 jx> CO I /\ I 2 ^ 

a = E(r (x,y)) = [t (v,(jj)| $ (v ; X,Y) dvdw 

d a — K 3 . 


(73) 


.1/2X .1/2Y 
-1/2X -1/2Y 


$^(V,w; X,Y) 


dvdco 


(74) 


where E denotes the statistical expectation operator. Due to the 
independence of the components of r(x,y). 


2 2 2 2 
a = E(r (x,y)) =0+0 
r s a 


( 75 ) 


Now, note that the effect of the detector is one of averaging 
or low-pass filtering (e.g., see (28), (45)). Although this is 
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necessary to preserve a practical signal-to-noise ratio, it 
produces some loss in resolution, or blurring. Similarly, the 
electrical filter which reduces the amount of additive aliasing 
noise also produces some blurring due to its averaging effect. 

Thus, there is a trade-off between blurring and additive noise 
(aliasing, detector and electronic). The selection of design 
parameters must, therefore, be made so that both the loss of 
resolution as well as the additive errors (accuracy) remain within 
acceptable limits. 

To find a quantitative means for judging the amount of blurring 
note that, in the frequency domain, the loss of resolution occurs 
when the higher frequency content of the signal is attenuated, which 
in turn reduces aliasing error. 

2 2 .00 ^ 

= E(M (x,y)) = 4)j^(v,w) dvdo3, 

= E(r^(x,y)) = |t^(v,w)t (v,oj) I ^ $^(v,oa) dvdoo. 

2 

Thus, cr^ is the variance that remains in the signal part of the 

2 

reconstruction, whereas is the variance of the original signal. 

Hence, the ratio a /o^, provides an indication of the amount of 

s M 

variance still remaining in r^(x,y). When the total frequency 

/N /V 

response T (v,Cjj)t(v,Cjo) is real, monotonically non-increasing, and 
no greater than unity in absolute value, this ratio provides a use- 
ful measure of blurring. A more general measure of blurring can be 
obtained as follows. Define blurring as the difference between the 


( 76 ) 

( 77 ) 
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original signal and the portion of the reconstruction due to the 


original signal, i.e., 

b(x,y) = r^(x,y) - M(x,y) 

Then, 

b(V,oj) = (t (v,(u) t(v,oj) - 1) M(v),w) 

K 

o 9 ooco,/^ i2^ 

= E(b (x,y)) = |Tj^(v,03) t(v,w) - l| 4 )^(v, 03) dvdw 

Note that, as defined in (78), blurring is the component that would 
still remain when the signal is sufficiently sampled, and no 
aliasing noise is present. 

r(v,oj) = M(v,o)) + (t_(v,oo) t(v,co) -1) M(v,w) + r (v,o); X,Y) 

ix 3 

r(x,y) = M(x,y) + b(x,y) + r (x,y) 

3l 

The ratio CJ, /a., provides a measure of the blurring, as a per- 
b M 

centage of the signal variance. The extent of aliasing error in 
the reconstruction can be compared by the ratio o^/o^. Finally, a 
quantity of significance is the signal- to-noise ratio a /a of the 

S d. 

signal variance to aliasing noise. 

It is important to note that both blurring and aliasing errors 

in the reconstruction depend on the particular method of recons- 

truction used, i.e., on the choice of T (v,oa) as can be seen from 

K 

(73), (77) and (80). Thus, the particular method of reconstruction 


(78) 

(79) 

(80) 

(81) 

(82) 
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can increase or reduce the errors involved, and should be selected 

with appropriate care to its effects on blurring and aliasing as 

well as other sources of error such as detector and electronic noise, 

uncertainty in knowledge of the sensor response parameters (the 

angular characteristics S(Xj^), the directional model R^(0', 

involved in t,(v,u3). No attempt was made, in this study, to 
d 

optimize the reconstruction method. The choice of the frequency 
response in (66) was made as it tends to neither reduce nor increase 
these errors, thus providing an indication of the extent of these 
errors in the sampled signal. However, to see the effect of 
further smoothing of the signal, the following reconstruction 
filters were also considered 


Trj^(v,03) = II(Nyv, Nyw) 


1, |Nyv| < Jg, |Nyw| < 

0, otherwise 


(83) 


The corresponding reconstruction r^(x,y) is given by 


r^(x,y) = ^ 

N Y j=-<» k=-' 


OO 00 

\ ^ \ / ^ x-iX . y-kY 

/ ^ ^ s(x^,y^) sin e — ■ - sincr 


NY 


Ny ’ 


(84) 


where Y is the effective diameter of the footprint. The variances 
of the smoothed reconstruction, r^^(x,y) , of the aliasing error in 
the smoothed reconstruction, and of the signal part of the smoothed 
reconstruction can be obtained in a similar way as for the unsmoothed 
case. Using (71), (73) and (75) we obtain 
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4 = = Lywy 


dvdw 


= ilS ^IS 


4=^4<->y» = 4 + 4- 


( 85 ) 

( 86 ) 
( 87 ) 


Three types of detector frequency responses were considered. 

These correspond to apertures of Gaussian, circular and diamond 
shapes. The patch considered was assumed to be at nadir, and the 
angular characteristics of the detector were assumed to be uniform 
over the IFGV. 

S(x,i|^) = 1 , cosx 1. (88) 

Thus, the point spread function of the detector was assumed to be 
unity within the footprint, and zero outside the footprint, except 
for the Gaussian case. In order to compare, the results from 
different types of detectors, the point spread functions were 
normalized as to satisfy 

lo. Lao dxdy = ^ (89) 

for each case considered, where Y is the diameter of the footprint 
for the circular aperture at nadir. 

2 2 2 

Y = h 6^ (90) 

where 6 is the IFOV, in radians. The spatial frequency responses 
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are given below. 


a) Gaussian: 
T^(V,0)) 


2 2 , 2 ^ 2 . 
-IT a (v +03 ) 
e 


, Si ^ y/2 


b) Circular: 

^ J (2TTap) 9 9 9 

T (v,o)) = — — , a = y/2, p = V + 03 

a Trap 

c) Diamond: 

Tj(v,03) = sinc(aV+bo3) sine (av-bo3) , a = , b = 

d 4 Z 


Perspective plots of these frequency responses are shown in Fig. 3. 
The electrical filter frequency response was selected as 


T (03) 
e 


1 ^ 11^1 

1 - (yo3) , |03| < - 

0 , otherwise 


Although this is not the frequency response of a causal electrical 
filter, it is representative of the magnitude characteristics of 
realizable filters. However, the phase characteristics that would 
introduce further blurring are not represented appropriately by 
this model. More realistic filters which operate in real time were 
modeled and coded in a simulation as rational transfer functions 
up to eight order, so that the capability of simulating the effects 
of these filters was developed. 

Finally, the radiant exitance at the top of the atmosphere was 
modeled as a random process (or field) which is homogeneous and 


(91) 


(92) 


(93) 


(94) 
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isotropic, i.e., wide-sense stationary. Since the particular 
values of the radiant exitance at the top of the atmosphere are 
not known a priori, and it is desirable to design the radiometer 
for a large class of possible signals, it is usual and pratical 
to model the radiant exitance as a random process. We will con- 
sider M(x,y) to represent the fluctuation of the radiant exitance 
about its mean value, say M, as the effect of a mean value M is 
to introduce a mean value in the output signal of the radiometer, 
as can be seen by 

OO 00 

Zoo /oo (M(x’,y’) + M) t(x-x', y-y’) dx’dy’ 

OO OO ^ 

= / f M(x’,y') T(x-x',y-y’) dx’dy' + M // T (x' ,y ’ )dx' dy ' (95) 

—00 —00 ^ QQ 


= s(x,y) + s 

(96) 

00 

"s = M // T (x’ ,y ' ) dx’dy’ 

(97) 


Since the effect of an average value can thus be taken into account, 
we assume that 

E(M(x,y)) = 0 (98) 

2 

(j)j^(x,y) = E(M(xQ+x,yQ+y) M(xQ,y^))= e (99) 

r = (xVy^)^. (100) 

The Wiener spectrum corresponding to the covariance function, 
cj)^(x,y), is given by 
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( 101 ) 


$fj(V,(D) 


2tt y 


2 

r 


O 


2 

M 


r , 2 2,2^ 2.13/2 

1 + Air y^(v +03 ) 


The Wiener spectrum in (101) can be obtained by assuming that 

M(x,y) is a random set of two-dimensional pulses whose width obeys 

a Poisson probability law with a mean width and whose magnitude 

2 

obeys a Gaussian probability law with zero mean and variance 
[15], [16], [17]. This interpretation provides a physical signifi- 
cance to the parameter It is important to note that by varying 

y , it is possible to distribute the frequency content of the signal 

with ease. As y decreases, the spectrum becomes flatter, favoring 
r 

2 

the higher end of the spectrum. When 0^^ is kept constant, the 
variance of the signal remains unchanged, so that meaningful com- 
parisons can be made. Finally, note that for the Lambertian case, 

L^(x, y, 6’, ({)’) = ^ M(x, y) , 

so that the analysis for the radiance would be essentially the 

/S 

same if the subscript M is replaced by L; i.e., by*^L * 
by $^(v,co) , etc. 


B. Information Density 


In situations where information in the input signal (e.g.. 


( 102 ) 
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the radiant exitance at the top of the atmosphere) is transmitted 
through a system (the radiometer) resulting in an output signal 
(the sampled digital radiometer output) , a criterion of interest 
is the density of the information transmitted. This criterion, 
which was introduced by Shannon [IB] for communication systems, 
was investigated for image quality considerations by Fellgett and 
Linf oot [19] and for aliasing effects by Huck and Park [20]. 

Since transmitting a signal which is known a priori does not 
increase our knowledge, the mathematical concept of information 
is based on the amount of uncertainty reduced by the transmission 
of a signal. The entropy [iS] of a signal is a measure of the 
uncertainty of receiving it. For a band-limited, second order 
stationary, Gaussian random process, g(x,y), with Wiener spectrum 
$ (v,cj), the entropy can be shown to be [l8, pp. 91 -95] , [20] , 

s 

H(g) = ly /2 -u /2 ?g(v,u)) dvdco. 


(103) 


where the spectrum ? (v,u) is limited to the band 

s 

{(v,oj)| |v| < Vq/2, |co| < (Dq/2}. 

In the previous section, it was shown that the reconstruction 
r(x,y) can be written as the sum of a signal part, r^(x,y), and an 
aliasing error part, r (x,y) (see (70)). The effects of detector 
and electronic noise as well as quantization error effects were 
neglected. Now, suppose that the detector and electronic noise 
from the amplifier and electrical filter, referenced to the input 
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of the filter, are modeled by Gaussian noise. Let n(x,y) be the 
noise added to the detector output when the scan point is (x,y). 

Thus, the input to the electronics can now be expressed by m(x,y)+n(x,y) . 
To simplify the formulation, we shall assume that the noise, n(x,y), is 
white when y is kept constant and x is allowed to vary. When the 
temporal power spectral density function of the noise is relatively 
flat, or alternately, when the temporal correlation function of the 
noise decays rapidly (i.e*> if the correlation function is almost zero 
in a fraction of the time required for one scan) , then the assumption 
is realistic. Since two scan points (x^,y) and (x^jy) will correspond 
to times t^ and respectively, separated by an interval, - t^, 

which, on average, equals the time required for one scan, n(x^,y) and 
nCx^jy) will be almost uncorrelated; hence, in the x direction the noise 
will appear to be white. Thus, let 

E(n(x+x' ,y+y ' ) n(x,y)) =4>^(y’) 6(x’) (104) 

The Wiener spectrum of the noise will then be dependent on one variable, 
as it is constant in V. 

/N 

4) (v,o)) =4) (co) . (105) 

n n 

As shown in the previous section (see (52) - (58)), the effect of 
the electrical filter is to process the noise only in the tD-axis. 

Adding the effects of sampling, A/D conversion and the reconstruction 
filter, the reconstructed signal can be expressed by 
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(106) 


r(x,y) = r (x,y) + r (x,y) + r (x,y) + r (x,y) , 

S a XI ^ 


where r^(x,y) and r^(x,y) are the parts of the reconstructed signal 
due to noise and quantization (A/D conversion) errors, respectively. 
The spectrum of the noise term can be found to be 


(f)^(v,a)) 

n 


uu 

E 

k=-oo 




n 


(O)-^) = |xv|<i |Yu|<i 


n 


(107) 


If the A/D conversion consists of rounding to the nearest allowable 
level, and is the maximum error that can occur (except for overflow) , 
then the quantization error can be modelled by 


$ (v,u) = E^/3. 

r 0 

q 


If the digital words consist of ri bits plus a sign bit, the 
number of allowable levels, k, is 


(108) 


K 


2 


ri+1 


Finally, if the signal is expected to have a dynamic range of j^-M, M J , 

(f>^(v,(D) = Eq = M/k. (109) 

q 3 k 

In the Appendix, it is shown that for a second order stationary 
real Gaussian random process, the aliasing error and signal processes 
are second order stationary and statistically independent. We further 
assume that both the quantization error and noise terms are Gaussian 
and statistically independent from the radiant exitance signal. The 
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Wiener spectrum of the reconstructed signal r(x,y) in (106), is then 
given by 

= $^(v,w) + $^(v, 0 j) + ?^(V, 0 J) + $^(v,co). (110) 

s a n q 

For the channel described by (106), where r^(x,y) is the input 
signal and r(x,y) is the output, due to the statistical independence 
of the additive terms, it can be shown that [l7, p. 99 ] the information 
density is given by 


h. = H(r) - H(r + r + r ). 

1 a n q 

Using (103), and manipulating 


h. 

1 


1 ^2X ,2Y 

2 - 

2X 2Y 


log2 



4>g(v,w) 


$ (v,oj;X,Y)+(i) (oj)+M^ic ^/3 

3 r 

n 


dvdo3 


(v, 0 )) = |t(v,(jo)| 




$^(v,co;X,Y) 



j=_oo k=-» 
(j ,k)ii^(0,0) 


( 111 ) 


( 112 ) 


(113) 


(llA) 


Note that when the cut-off frequency of the low-pass electrical 
filter is chosen so as to attenuate the frequency content outside the 
Nyquist band significantly, the aliased portion of the noise becomes 
small relative to the portion of noise within the passband; i.e.. 
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E 

k=-“ 

k^O 


Te(U3-|) 


<1,^ («)--) « 


T (w) 
e 


(p (w), 

n 


Yw 


(115) 


2 

(116) 

n 

Thus, the information density transmitted through the system can be 
computed using (112) - (116). To see the effect of specific parameters 
more explicitly, consider the following normalizations: 




(117) 


t(v,w) = K t'(v,co), t (co) = K t’(o)), (118) 

T e e e 


where the prime " ’ " denotes a normalized quantity, is the total 

"d-c" gain of the radiometer system, the d-c gain of the electronics, 

2 2 . 

and the variances of the input signal and the noise, respectively. 
Substituting these expressions into (112), 


, 1 .2X .2Y T 

i " 2 -i -A. 

2X 2Y 


1 + 


t'(v,w)| $j^(V,0)) 


(v,to;X,y)+ |t'(uj) 1 $ (u)+ 

a \o/\K/ e n 

\ n/ \ e/ 


-2 

M K 
2 2 
M T 


dVdw 


(119) 


The expression in (119) shows the dependence of the information 
density on various system parameters explicitly. For example, if the 
total system gain, K^, is increased by increasing the electronic 
amplification K^(i.e., remains unchanged), the information density 

would remain unchanged, assuming that the dynamic range of the signal 
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I 


increases proportionately. On the other hand, if is increased by 
increasing the gain of the detector or the lens configuration, the 
effect of the noise would be reduced, while the aliasing error term 
would remain unchanged, assuming remains unchanged. To see this more 
clearly, let t>e the signal to noise ratio of the input signal 

(radiant exitance) to the noise (detector plus electrical) both referenced 
at the output of the detector; then 


K 

s/n K M n 
e 

On the other hand, the standard deviation of the signal at the 
input of the A/D converter is 


( 120 ) 



K‘ 


+ a 


n 



cj K if a / >> 1. 

M t’ s/n 


( 121 ) 


Since the dynamic range of this signal would be directly proportional 
to its standard deviation, let 


M = c K 
M T 

Substituting (120) and (122) into (119), 


( 122 ) 


h. 

1 


A. 

1 .2X ^2Y 

2 - 

2X 2Y 


log2 


1 + 


t’(v, 0))| (|)jJj(v,a3) 


dvdw (123) 


(v,to;X,Y)+ a / |t' (U 3) V (t^)+ — K ^ 

a s/n' e ' n 3 


It is seen that the combined effects of the gains and 0„, a can be 

M n 

combined into the parameter i.e., the signal to noise ratio 

referenced at the detector output. The effects of shaping the transfer 
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functions t'(v,w) and f ^ (w) are less explicit and have been investigated 
in the previous section for aliasing and blurring. The effect of the 
sampling intervals is shown in Figure 10. 
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III. RESULTS 


The models described in section II were used to make a para- 
metric study of the effects of various design parameters such as 
aperture shape, continuous versus step scan, sampling interval, 
etc. , on aliasing error. A computer program which computes the 
variance of the signal and aliasing error, and signal-to-noise 
ratio was developed. The program generates the Wiener spectrum 
of the optical signal, the frequency response of the lens, the 
frequency response of the selected detector aperture shape, and 
the transfer function of the electrical filter. Then, using these 
functions, the Wiener spectra of the continuous and sampled radio- 
meter output are computed. The spectrum of the sampled signal was 
generated by including the effects of the eight sidebands closest 
to the fundamental band at the origin, as these have the largest 
effect on the spectrum, while further sidebands have negligible 
effect within the reconstruction band {(v,03)[ ]xv| < h, |Yo)[ < ' 2 } . 
In all the runs made the sampling intervals in the along- track and 
scan directions, x and y respectively, were chosen to be equal; 
i.e., X = Y. 

To obtain results which are useful in a more general context, 
the independent variables x and y were normalized by the diameter, 
Y, of the effective circular footprint. Thus, the values X/y=Y/y= 1 
correspond to the contiguous sampling case for a circular IFOV, 
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irrespective of the particular value of Y* Similarly, the average 
pulse width of the signal is also normalized, so that the normalized 
parameter U^/y represents the average number of zero-crossings 
occurring within the footprint. To make the results applicable to 
different levels of signal power, the variances of the signal and 
aliasing error were computed as ratios of the variance of the in- 
coming radiant exitance. As the system is linear, any particular 
signal level can be readily obtained from the ratios. The integrals 
required to compute the variances of the signal and aliasing error 
(e.g., (72), (74)) were evaluated using simple numerical integration 
techniques . 

Figure 6 shows the normalized variances of the signal and 
aliasing error components of the reconstruction given in (65) and 
(70) and the signal- to-noise ratio as a function of the normalized 
sampling interval X/y, various Wiener spectra as input radiant 
exitance; the spectra are parametrized by y^/Y* The three cases 
shown correspond to the Gaussian aperture with step-scan, the 
circular aperture with step-scan, and the diamond aperture with 
continuous scan. 

Figure 7 shows the normalized variance of the aliasing error 
in the smoothed reconstruction given in (84) and the signal-to- 
noise ratio as a function of the normalized sampling interval for 
varying amount of smoothing indexed by N. 

Figure 8 shows the aliasing error variance and the signal-to- 
noise ratio for the smoothed reconstruction when continuous scan. 
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hence electrical filtering, is used. To see the effect of continuous 
versus step-scan to aliasing error, compare with Figure 7. 

Figure 9 shows the normalized variances of the signal and 
aliasing error parts of the smoothed reconstructions, and the 
signal-to-noise ratio as a function of the frequency content of 
the input radiant exitance indexed by y^/y. It is seen that while 
smoothing reduces aliasing error, it also increases blurring as 

measured by a /o . 

S i-j 

Figure 10 shows the information density transmitted through 
the radiometer for the unsmoothed reconstruction as a function of 
the normalized sampling interval. Plots for a radiant exitance 
Wiener spectrum corresponding to white noise with unit power, as 
well as various levels of coloring are given. Note that the white 
noise spectrum has infinite variance, so that the results are not 
readily comparable. 
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IV. CONCLUSIONS 


1. The numerical results substantiate the observation by 
Mertz and Grey [2j and Schade |]3j that all reasonable spot 
intensity profiles and photosensor aperture shapes of equivalent 
size result in about equal blurring, but that some profiles and 
shapes tend to suppress aliasing better than others. The Gaussian 
spot intensity profile of television cameras is appreciably 
superior to the circular photosensor aperture commonly used in 
optical-mechanical scanners. However, aliasing with optical- 
mechanical scanners can be effectively suppressed by careful 
photosensor aperture and electronic filter response shaping. 

2. Aliasing errors can be the most significant source of 
degradation in contiguous reconstructions of a discrete signal. 
Consider, for example, a discrete signal from an optical-mechanical 
scanner with a circular photosensor aperture and contiguous coverage 
(i.e., X/y = Y/y = 1), and a random radiance fluctuation with a 
mean spatial detail equal to the sampling interval (i.e., y^/y' = 1) . 
The resulting rms signal-to-noise ratio a /a is then only 5, 

S 3. 

whereas imaging systems and radiometers are commonly designed to 
provide signal-to-noise ratios of two to three orders of magnitude 
in typical applications. 

3. Aliasing errors tend to decrease rapidly with decreasing 
sampling intervals. Improvements of one or two orders of magnitude , 
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depending on the spot intensity profile or aperture shape, can be 
attained by decreasing the sampling intervals X/y = Y/y from 1 to 
0.7. Still further decreases in sampling intervals continue to 
provide rapid improvements for the Gaussian profile, but relatively 
gradual^ improvements for the circular and diamond apertures . This 
difference in performance between the Gaussian profile and the two 
apertures is directly attributable to their frequency response 
characteristics: the frequency response of the Gaussian profile 

decreases monotonically with frequency, whereas the frequency 
response of the two apertures exhibit sidelobes. The reason for 
the improved performance of the diamond over the circular aperture 
is that the sidelobes of the diamond aperture are suppressed along 
the two sampling directions (i.e., by its shape along the co- 
direction and by an electronic filter along the V- direction) . 

4. The improvements that can be attained in suppressing 
aliasing by photosensor aperture shaping and decreasing the sampling 
interval tend to be independent of the statistical properties of 
the radiance fluctuation. 

5. Spatially smoothed reconstructions reduce aliasing error 
at the cost of increased blurring. But there is a net loss in 
signal-to-noise ratio (i.e., O /o O /o ), and hence loss of 
information about spatial detail. This suggests that the pass- 
band of the reconstruction filter should generally be equal to 
the sampling passband. Narrower reconstruction filters (i.e.. 
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smoothed reconstructions) should be used only if spatial detail 
is of secondary importance to the accurate estimation of absolute 
radiance magnitudes of coarse detail (as it is commonly in the 
reconstruction of radiometric measurements) . 
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APPENDIX A 


SAMPLING OF SECOND ORDER STATIONARY RANDOM PROCESSES 


Consider a second order random process {x(t), — «» < t < , 

defined on a fixed probability space A, P) . Let 

E(x(t)) = 0 

E(x(t)x*(t')) = 

Suppose that the random process is sampled at intervals T, 
to obtain the sequence {x(kT) , - “ < k < °°} . Whatever the physical 
nature of the sampled process (whether a sequence of numbers stored 
in a computer, the output of a sample— hold device, etc.)> it con- 
tains information about the continuous process x(t) . In many 
practical applications , it is desirable to obtain a continuous 
random process, x^(t) , from the samples of another random process, 
x(t) , such that > resembles x(t) is some fashion; e. g., 

E|x(t) - x^(t) I is small or minimum. Such a process, x^(t) , will 
be called a reconstruction of x(t) . In the following, we shall 
assume that the random processes considered have zero mean with- 
out loss of generality, as the mean (if finite) can be subtracted. 

Def. 1 

A random process is said to be second order stationary (or 
wide sense stationary) if it is of second order and its covariance 


1 


* denotes the complex conjugate. 



depends only on the difference of its arguments, i. e.. 


E(x(t4T)x*(t)) = t, TeR. 


Second order stationary processes have harmonic (or spectral) 
decompositions as stated in the following theorem (for a proof see 
Loeve , pp. 474 - 483). 


Theorem 1 

A second order random process x(t) with tcR has a harmonic 
orthogonal decomposition 


x(t) = e^^'^dX(w) 


where {X(co) , (jOeR) is a process with orthogonal increments such 
that whenever 


E\X(i^P - X(w^) 1^ = - F(u3^), 


and F(w) is of bounded variation if, and only if, 

x(t) is second order stationary and continuous in quadratic mean 

(q.m.) at one point. Then, 


li'S' il x(t)e = x(h)f) - X(u-) 


l.i.m. 1 . e 

T 0° 27T -T 


-io!t -i(o3+h)t 


dt = X(cjofh) - X(co) , 


^ 1 
2 


X(w) = 4 X(ojf) + X(oo- 


R denotes the real numbers 
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The theorem remains valid if the argument t is restricted to the 
integers, the continuity of x(t) is deleted, and the integrals in 
(6), (7), are replaced by appropriate sums. 

The stochastic integrals above should be interpreted as 
Riemann-Stieltjes integrals in q.m. (see Loeve |2lJ , p. 472), and 
l.i.m. denotes the limit in quadratic mean. 

For a second order stationary process, x(t) , continuous in 
q.m., we shall call any process X(o}) satisfying (4) an Integrated 
Fourier Transform (IFT) of x(t). A non-decreasing function F (oo) 
satisfying (5) will be called an Integrated Power Spectrum (IPS) 
of x(t). If F(o 3) is absolutely continuous, then the derivative 

S(u) = -^(u) (9) 

3 

is called the Power Spectral Density (PSD) of the process x(t) . 

As in the case of deterministic signals, the concepts of power 
spectra (or harmonic decompositions) and Fourier transforms play 
an important role in the analysis of random processes. 

Def. 2 

Let (x(t) , teR} be a second order random process, and 
{t|,, k = . . ., -1, 0, 1, . . . } an increasing sequence of sampling 


More generalized definitions applicable to non-s tationary process 
are possible and will be studied later. 
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points. Then {x^(t) , teR} is a real (complex) linear reconstruction 
of x(t) if there are real (complex) functions such that 

n “ 

x^(t) = 

In particular, if tj^ = kT, the linear reconstruction 

OO 

x^(t) = ^S^x(kT) sinc(k - t£R (11) 


will be called the Shannon reconstruction of x(t) if the limit in 
(11) exists. 

Note that, for the Shannon reconstruction. 



jl, k = j 
|o, k j ’ 


( 12 ) 


X (t, ) = x(t ) , t = kT. 

s k k k 


(13) 


Condition (12), which implies (13), is clearly a desirable one 
when trying to reduce the difference between x(t) and the recon- 
struction. 


Def . 3 

Let {x(t), t£R} be second order stationary and continuous in 

til 

q.m. , with IFT X(o)) . Then the k sideband of the process for a 
uniform sampling interval T is 


Xj^(a3) = X(aj+ka)Q) - X(kooQ), = 27 t/T, k = 0, ±1, . . . (14) 

th 

The k sideband signal of the process is 

^ /2 

k = 0, ±1 (15) 
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Lemma 1 


Let {x(t), teR} be second order stationary and continuous in 
q.m. Then the sideband signals {x^(t) , t£R} are uncorrelated, 
second order stationary processes; the sidebands 

{Xj^(w) , ~^q/2 < w < ^J^q/ 2} are processes with orthogonal increments 

such that 


E([X|^(ft)2) - X|^(to^)][Xj(u3^) - Xj((Oj)] *) = 0, k i* j (16) 

whenever -w /2 < co , w , co 03 < o) /2; and 

U 1 2 3 4 0 

E(|x^(0)2) - 1^) = ^1 - ^2 

F^(oj) = F(03 + k03^) . (18) 


Proof : 

First note that (17) follows from (14) and (5); since, whenever 
^ “2 

1^) = E(|x(u^+ku^) - X(w^+ku^)|^) (19) 

= F(032+1<O3q) - F(03^+k03Q). (20) 

Now, let 03^£[-03 q/2, o)q/ 2), i = 1, 2, 3, 4; then 

E([x(o32+ko3Q) - X(o)^+ko3Q)] [x(o)^+jojQ) - X(oj^+jo)Q)] ) = 0, k j (21) 

since X(03) has orthogonal increments (Th. 1) and the intervals 
[min(O3^,0J2)+k'03^, max ( , CO 2 ) +k 03^ ) and [min(o3^,03^)+ja3Q, max(oj^ ,03^)+jo3Q) 
are disjoint if k j . Substituting (14) into (21), (16) follows. 
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To show that the sidebands have orthogonal increments, suppose 


/2 < w w < ^./2; then 

0 1 2 3 4 0 


E([\(w2) C V"^4^ *)=E([x(o)2+hi^Q)-X ^^^+^0)^)] 

[X ( (jOq ) - X ( oj^'^kioQ ) ] ) =0 

Thus, the assertions about the sidebands have been proved. To see 
that the sideband signals are stationary, it suffices to note that 


w /2 . w /2 , 

R^(t,t')=E(x^(t)x*(t'))=E(/^^^2 

-\/2 ^ • 

Thus, the covariance R^(t,t') depends only on the difference of the 
arguments, i.e., t-t'. 

? w /2 

E(|xj^(t) I )=C/2 dF^(“) = \(“o/2)-Ek(-“o^2) < «>, tcR. 

Hence, is second order stationary. 

Now let k ^ j and t, t' eR, 


E(x^(t)x*(f)) = ^V“> 4"/2 


iojt 


/ 2 • 1 ? f, 

0 dx*( 0 i*)) 


= 0 


where (27) follows from (16) and the definition of the integral. 
Hence, the sideband signals are uncorrelated, and the lemma is 
proved , 


( 22 ) 

(23) 

(24) 

(25) 

(26) 
(27) 
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Def . 4 


Let (x(t) , teR} be second order stationary and continuous in 

ttl 

q.m. with IFT X(w) . Then the k two-sided sideband is 

X^(o3) = X^(co) + k = 0, ±1, ±2, . . . . 

th 

The k two-sided sideband signal is 

3^(t) = x^(t) + x_^(t), k = 0, ±1, ±2, . . . . 

Lemma 2 

Let {x(t), teR} be second order stationary and continuous in 
q.m. with IFT X(w). If a < b and the complex function, f, is 
continuous on [a, b], then 

f (oo)dX(o))=/^ f (oa)dX(w)+f (b) [x(b)-X(b)]-f (a)[x(a)-X(a)] 

S. 3 . 

where X(w) is defined by (8). 

Proof : 

Since X(w) coincides with X(w) when CO is a continuity point 
(in q.m.), the second order process 

X(oj) = X(w) - X(co) 

vanishes except for a countable number of discontinuity points. 
Hence, X((j 0) is zero on a set C dense in [a, b) . Now choose a 
sequence of finite partitions, P^, of [a, b) , whose points lie in 
C with the possible exception of a and b. Since C is dense, the 
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( 28 ) 


(29) 


(30) 


(31) 


sequence can be chosen so that the norm of the partitions converges 
raonotonically to zero. Then, 

f^(a))dX(oa)=F(aj )X(b)-f(w. )X(a), (32) 

a P nn In 

n 

where is in the first subinterval of P and w is in the last, 

In , n nn 

and fp(ci3) is the restriction of the continuous and bounded function 
n 

f((jo) to P . Since ca. a, w -> b and f is continuous, 
n In nn 

f (o))dX(aj)=f (b-)X(b)-f (a+)X(a) . (33) 

a 

Hence, the integral on the LHS of (33) exists; since f(cij)dX(w) 

0 . 

exists, (see Loeve [21j , p. 473), then so does / f(oj)dX(ca); and 
the assertion is proved. 

In particular, note that if X(o)) is continuous in q.m. at a and 
b, then the last two terms in (30) vanish and the two integrals are 
equal |a, ej. It is also interesting to note that for a continuous 
integrand, if w is an interior point of discontinuity, then the 
integral depends on the jump X(a}+)-X(co-) , but does not depend on 
X(to), as evidenced by the preceding lemma. 

Lemma 3 

Let {x(t) , teR} be second order stationary and continuous in 
q.m. with IFT X(w) . Then, the second order process X(oj) , defined by 

X(co— ) , U) > 0 

X(oo) = (34) 

X (Ci3+) , (jO < 0 
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has orthogonal increments; furthermore, 


x(t) = / e^^^dX(o)) , teR 


( 35 ) 


Proof : 

Consider four points on the real line such that 
-oo < < 0 )^ < co^ < co^ is necessary to show that 

__ ^ 

E([X(w^)-XCa)3)][x(aj2)-X(a3^)] ) = 0 (36) 

First let > 0, then 

X(oo )-X(go )=X(oj,-)-X(to -)= l.i.m. X(o) -h)-X(oo -h) (37) 

^ ■^h4'0 ^ 

Now, if (^>2 - then 

lim E([x(oo -h)-X(oo -h)] [X( 0 ) -h)-X(o) ±h)]*) = 0, (38) 

hiO ^ ^ 

if 0 O 2 < 0, then 

lim E([x((jo,-h)-X(co -h)] [X(u) +h)-X(o) +h)] ) = 0, (39) 

h'^O ^ z 1 

Hence, (36) holds when 00 ^ > 0. Finally, if < 0, 

lim E([x(w ±h)-X((i3 +h)] [x(co +h)-X(w +h)]*) = 0. (40) 

h'l'O z 1 

Hence, X(o3) has orthogonal increments. 

To show (35), recall that F(u)) is of bounded variation; hence, 

the following limits exist. 

l.i.m. X(±oj) = X(±o°), (41) 

CJ -> CX3 
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(43) 


0 < E X(-ta)-X(^£ji>f) = F(-Wf)-F(-w) ^ 0 


So that X(±°o) can be defined as 


X(±«o) = l.i.m. X(±w) = X(±“). 

0 > 00 


(44) 


Since X(w) and X(w) differ only on the countable set of discontinuity 
points, using Lemma 2, 


/.b icot V ^b iuit,^. V ibt~,, . , iat~. ^ 

/ e dX(u3)=/ e dX(w)-e X(b)+e X(a), 

Si 3, 


(45) 


X(0)) = X(w)-X(w). 


(46) 


From (42) , (43) we see that 


l.i.m. e^^^X(b)-e^^^X(a) 

a -oo 
b oo 


= 0 


(47) 


Thus, applying the limit as a and b ^ “ to (45), (35) follows. 


It is important to note that X(o)) has orthogonal increments if, 
and only if, X(cj) is continuous in q.m. on the real line, as can be 
seen by 


E([x(M^)-X(Uj)][x(l^)-X(tt^)]*) = 


0 


Oi) c(jl) 


^[f(lo.+)-F(o) -)], a) =0) 


( 48 ) 


2 3 


whenever < O)^. Thus, X(co) cannot be used in place of 



X(o)) , although it has symmetry about its discontinuity points. 


Lemma 4 

Let {x(t) , tcR} be a real, second order stationary process, 
continuous in q.m.. Let X(oj) be given by (34). Then, 

[x(a32)-X(o)^)]* = X(-03^)-X(-W2), ^ 0 

X*(w) +X(-cu) = Z, (JO ^ 0 

where Z is a real random variable. Furthermore, (jO is a continuity 
point if, and only if -(jo is a continuity point. 


Proof : 

Using (7) , we have 


ioJ,t i(jO„t 


r- - T* It 

X(to_)-X((jO.) =l.i.m. r X(t) ^ - dt 

L Z TJ 27T-T -It 


X OO 


i(jO t ico t 

1 - 1 / \ e -e ^ , 

=1.1. m. J x(t) : dt 

27T -T it 

X 00 


=X(-oj^)-X(-(i32) , co^oR. 


Since XCcxi) and X((iO) are identical except on a countable set 
of discontinuity points, (49) holds except on the discontinuities. 
As X(co) is continuous from the left for (jo > 0, and from the right 
for CO < 0, (49) holds by applying the appropriate limits. Rearran- 
ging (49), 


(49) 

(50) 


(51) 

(52) 

(53) 
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X*(a32) +X(-W2^=^* V ^1* “2 ^ ° 

so that Z does not depend on the value 03. To see that Z must be 
real, it suffices to set “1 “ -“2 

Now let ^2 = w + h, 00 ^ = O); then 

E|x(co+h)-X(o)) I ^=E|x(- 0 )-h)-X(-CA)) I oj 5 ^ 0 , (jD+h?^0. 

Hence, as h 0 (o)+-h^0) , either both limits converge to zero and 
X(co) is continuous in q.m. at W and -w, or (exclusive) X(^) is 
not continuous at either point. Thus, the assertion is proved. 

Theorem 2 

Let {x(t), teR) be a real, Gaussian, second order stationary 
random process, continuous in q.m.. Then the two-sided sideband 
signals {x^^(t) +x_j^(t) , teR} , k = 0, 1, 2, . . ., are statistically 
independent, second order stationary, real Gaussian random processes 
where 

X^(t) = /„°/2 s’- "^dXi^Cto), k= 0 , ±1 “0 ^ 

3t|^(u) = X(t«-ki.)Q)-X( k to ) toe[-to^/2, "q/2). 

The two-sided sidebands {Xj^(to) +X_^(to) , co£[-to /2, to /2)}, k=0, 1, . . 


(54) 

(55) 

(56) 


(57) 

(58) 


are statistically independent, Gaussian random processes with 
orthogonal increments. 



Proof : 


From Lemma 1, recall that the sideband signals are uncorrelated 
and second order stationary. Thus, 

E([xj^(t) [x (t') +x_ (t') ] *) = [Rj^(t-t' k, j=0,l, . . . , 

(59) 

So that the two-sided sideband signals are also uncorrelated and 
second order stationary. Hence, to obtain independence, it suffices 
to show that the processes are strongly Gaussian. First note that 
X(oj) and, hence, X^(co ) , are Gaussian; e.g. see Loeve, pp . 485-486. 


x^(t) = 


0^0/2 . 

■ ol/2 . 


* 00/2 

,0' -ioot-* 

-w /2 ® 

0 


dX^(U)) 


■ V'" 


_ w /2 

O' _-i£0t— , ,0 lwt,T7 




(60) 

( 61 ) 


where we have used Lemma 4, and the continuity of e 


iwt 


To see this 


in more detail, consider the partition ~co^/2 = w < to < ... <to = (i) /2. 

012 n+1 0 

The Riemann-Stielt jes sums are 




(00 +k(^_)-X(03.+kW_) 

j+1 0 3 0 J 


(62) 




i (- w’. ) t 


[^X(-a3^-k<^Q)-X(-w^_^^-ka)Q)j (63) 




j=l 

^ ifjj't, 






(64) 

(65) 
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where {oj , j = 1, 2, . . n+l} is the partition 

j 

-o)/2 = w = -w , . < w- = -w < ... < C = - 0 ) = to /2 

0 1 n+l 2 n n+l 0 


^n+l-j’ ^ 2, . 


. . , n. 


As the norm of the partition approaches zero, the sums on the LHS 
of (62) approach x^t) , while (65) approaches x_^^(t) . Thus, 
x^(t) + X ^(t) is real, hence strongly Gaussian, and independence 


follows. 

From Lemma 1, recall that the sidebands are uncorrelated and 
have orthogonal increments, hence the sum of any two and, in 
particular, the two-sided sidebands have orthogonal increments. 
Now, let 


= Xj^(u)) + -u^/2 < u < U)p/2. 

Then, 

Y*(w) = 

E(Yju) Y.(l»')) = -E(Y (10) Y*(-“')) = 0, k # j 

Thus, the two-sided sidebands (Y^(to) , aje|^-0jQ/2, 0JQ/2)}are 
strongly Qaussian and uncorrelated; it follows that they are inde- 
pendent, and the proof is complete. 

It is of interest to note that {X(w) , coeR} does not have 


( 66 ) 

(67) 


( 68 ) 

(69) 

(70) 


independent increments except for the trivial case where it is constant 
except for a possible jump at CjO=0, which corresponds, to x(t) having 
sample functions which are constant. This can be seen from 



(71) 


E([^X(W2)-X(ti3^)] [x(-w^)-X(-'032)]) = F(W2)-F(o3^), 

whenever 0 < The following lemma will be needed for 

Theorem 3. 


Lemma 5 

Let {X (oj), a<w^b}n=l, 2, . . .be processes with 
n 

jointly orthogonal increments such that 


e|x (co.) 

n 2 


- X (W ) 
n 1 


= F^Co.^) - F^Ccop , 0,^ 


^“2 


X(o)) = l.i.m. X (to) , a < (0 < b 

n ^ oo ^ 

Let f be a bounded complex function on [a, b] . If the integrals 
below exist, then 

f(co) dX(u)) = l.i.m. f(oji) dX (to), 
a ^ oo a n 

and {X((jo), a < co < b} has orthogonal increments. 


Proof : 

Let A.X = X(oj.) - X(o)!), i = 1, 2. 
in n 1 n 1 

Then, whenever to < to^^ £ (jc^ < , 

E(A,X A,X*) = lim E(A„X A-X*) = 0. 

2 1 ^ ^ ^ 

Hence, X(co) has orthogonal increments. Since X^(io) , have 

jointly orthogonal increments, i.e.. 


(72) 

(73) 

(74) 

(75) 

(76) 
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( 77 ) 


E(A,X„ AX*) = 0, 
Z Tn in 

it follows that 


E(A^X A,X*) = lim E(A.X A,X*) = 0. (78) 

2 In m-x» 2 m In 

Thus, by direct computation, it can be seen that the process 

X (oj) = X((jo) - X (w) 0, a < 0} < b, (79) 

n n 

also has orthogonal increments. Now, 

e|/^ f(w) dX(w) - f(oj) dX(to) 1^ = e|/^ f(w) dX(03) 1^ (80) 

a a n a n 

= (to) 1^ dF (03) (81) 

a n 

< dF ((jj) = M^[f (b) -F (a) ^ 0, (82) 

- a n n n 


where M is an upper bound on f, and 

F (U) ) - F (u ) = e|x (to.) -X (u,) 1^, to, S to, (83) 

n Z n 1 ' n 2 nl' l 2 

Hence, (74) follows. 

Theorem 3 

Let {x(t), teR} be second order stationary and continuous in 
q.m. Then, 

n 

X (t) = l.i.m. S x(t) = ^n^t) + x (t),. teR, (84) 

r k=-m ^ ° 
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( 85 ) 


X 


(u) = 1.1. m. X,(IU) = Xlw) + X (U) , < CO < 

m , n-><» k=— Til K- u a z z 


w /2 

X (t) = fJ^/2 dX (tD), teR, 

0 


( 86 ) 


oo 

E|x^(ai2)-X^(u^) 1^ =.5„o F(‘»2+kUg)-r(“j^+kWQ) = ’ “l ^ 


E(x^(t) xJCt')) = / °/2 e^“^"‘"'^(jr^(u)). (88) 

X (t) (X (w)) and x (t) (X(co)) are uncorrelated. Furthermore, if 

(J vJ d. 3. 

x(t) is real and Gaussian, then x^Ct) and x^(t) are statistically 

independent, real, Gaussian processes; X^(to) and X (co) are statistically 

U 3 

independent, Gaussian processes. The processes x (t) (X(to)) will 

3 3 

be called the aliasing error signal (aliasing error IFT) . 


Proof : 

To show the convergence in (84) , recall that (Lemma 1) , the 
sideband signals are uncorrelated. Using (25), we obtain 

e|X) 1^ = S E|X|^(t) 1^ = 22 E, 

k=-m k=-m k=-m 



(89) 



-A ^ U. 

2 ^^0 





< F(o°) 


F(-oo) < oo,(90) 


which is sufficient for the desired convergence (e.g., Loeve [2lJ , p. 456). 

Similarly, for wef-o) /2, co /2), 

U 0 
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k=-m 


Ej \(w)| = e| ^ XCcjiH-koo^) - XCkOJ^) 

’ k^m 

n 

= ^ E|x(w+ka)^) - X(kw^) 


k=-m 
n 


= ^ iFCio+kUp) - F(kUg) 


( 91 ) 


(92) 


(93) 


n 


k=-m 


to. 


"o 


< 2 ^("T " Fi- F(“) - F(-“).(94) 


2 Q 


Thus, the limit in (85) exists, and X^(ui) is defined. 

9 — _ 9 

e|x^(o^2) - I = E V“2^ ■ I 

k=-00 

oo 

= E e|\(“2) - w I' 

k=-a> 

OO 

= 2 F(o32+k{jjQ) - F(o 3 +kwQ) 

k=-oo 

results in (87), where Lemma 1 has been used. Now, to show (86), 
let 


(95) 


(96) 


(97) 


X (w) 
rn 

n 

= E 

X^(co) 

X^(io) , 

(98) 

k=-n 




n 





II 

M 

x^(t) 

iLli.' x^(t). 

(99) 

k=-n 




n 

= E 

k”-n 

\/2 


(100) 


Using (16) in Lemma 1, note that 
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( 101 ) 


n 


m 


E(A_X A,X ) = E( Z) A-X. Z A^X*) 
2 rn 1 rm j' 


n m 


— * 


= £n ECA^X^ A^X.) 


= kfcp ""A) = 

whenever the increments (defined by 75) correspond to disjoint 

intervals, and p = min(n, m) . Hence the processes, X (w) , have 

rn 

jointly orthogonal increments, and Lemma 5 applies. Taking the 
limit as n of both sides of (100) results in (86). By Lemma 5, 
{X^(co), - — j < oj < ~~ 2 ^ orthogonal increments; thus. 


E(x^(t) x^t’)) 


= E( 


' ^ l2 

itil /2 ® 

o' 


■ oj /2 

e"^"'dX (c) 

0 


"I * 


- ioo(t-t’) .p /,.s 

- U^/2 ^ 


If x(t) is real and Gaussian, to show the desired properties 

it suffices to note that x^(t) , x (t) , X„(to) and X (oo) can all be 

U a 0 a 

expressed in terms of the appropriate two-sided processes and apply 
Theorem 2. 


x^(t) = igjxQ(t) +x_^(t)J , Xq(iji)) = l5j^XQ(t^J) +X_ q(oj)J . 

“kS • 

Since x (t) (X (oj)) and x (t) (X (uj)) are Borel functions of inde- 
U U a a 

pendent processes, they are themselves independent, and the proof is 


(102) 

(103) 


(104) 

(105) 

(106) 
(107) 
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complete. 


Corollary 1 

Let {x(t), teR} be second order stationary and continuous in 
q.m. If the IPS F(u}) is absolutely continuous on R, with PSD S(o3) 
then the PSD of the process x^(t) , say is given by 

00 

5(ciO) = Y2 ), -w /2 < CO <o) /2, 

^ k=-o® u u u 


Proof : 

From (87) we have 


n 


F(co„) -F(co^) = lim ^ / S (co+kw^) dco 

^ r -L n-K» ^ 


00 , 


= lim / 


n-x» 


00 , 


S(o>fkco ) 
k=-n ^ 


dco 


Since F(co) is non-decreasing, the PSD 5 (oo) > 0, oo£R. Hence, the 
integrand in (110) is non-decreasing as n t «5. By the Monotone 
Convergence Theorem, 


00 , 


= I 


£ S(U)+k03 ) 

k=-CX5 


dtJj; 


so that F^(oo) is absolutely continuous and its derivative is given 


by (108). 


(108) 

(109) 

( 110 ) 

( 111 ) 
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Figure 1 Satellite, Orbit and Top of the Atmosphere Geometry 








(c) DIAMOND 

Figure 3 Frequency Response of Detectors 
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(a) Gausjion spol inlonjify profile (without 
electronic tiller). 


(b) Circular photosensor aperture shape (without 
electronic filter}. 


(c) Diamond photosensor operture t'hape (with 
electronic filler). 


Figure 6 Parametric Aliasing Error Curves. 
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Figure 9 Aliasing Error Versus y /y for Smoothed Signals 
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